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C^") ' A deformed relativistic Hartree-Bogoliubov (DRHB) model is developed aiming at a 

proper description of exotic nuclei, particularly deformed ones with large spatial exten- 
sion. In order to give an adequate description of both the contribution of the continuum 
and the large spatial distribution in exotic nuclei, the DRHB equations are solved in a 
Woods-Saxon basis in which the radial wave functions have proper asymptotic behaviors 
at large distance from the nuclear center which is crucial for the formation of halo. The 
formalism and the numerical procedure of the DRHB model in a Woods-Saxon basis are 
briefly presented. 

•i-H . 

' Keywords: Relativistic mean field model; Bogoliubov transformation; Deformed halo; 

<-H ' Woods-Saxon basis. 

1. Introduction 

One of the exotic phenomena observed in nuclei close to drip lines is the halo in 
which the extremely weakly binding property leads to many new features, e.g., the 
coupling between the bound states and the continuum due to pairing correlations 
and the large spacial density distribution. In order to give an adequate theoretical 
description of the halo phenomenon, the asymptotic behavior of nuclear densities 
at large r must be considered properly and the discrete bound states, the contin- 
uum and the coupling between them must be treated self consistently. This could 
be achieved by solving the non- relativistic Hartree-Fock-Bogoliubov (HFB) 1-3 or 
relativistic Hartree-Bogoliubov (RHB) 4 ~ 6 equations in coordinate space which can 
fully take into account the mean-field effects of the coupling to the continuum. 
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In Refs., 2-6 the spherical symmetry is assumed. Since most of the known nuclei 
are deformed, whether or not there exist deformed halos and what new features are 
expected in deformed exotic nuclei are very interesting questions 7-12 which could 
be answered by the deformed counterparts of the HFB or RHB models. Neverthe- 
less for deformed nuclei, to solve the HFB or RHB equations in coordinate space 
becomes much more sophisticated and numerically very time consuming. Although 
many efforts have been made to develop non-relativistic HFB models either in (dis- 
cretized) coordinate space or in a basis with improved asymptotic behavior, 13-20 the 
deformed relativistic Hartree-Bogoliubov model has been developed only in conven- 
tional harmonic oscillator basis. 21-24 

The deformed relativistic Hartree equations have been solved in a basis in which 
the basis wave functions are calculated from solving the Dirac equation with spher- 
ical Hartree potentials. 25 ' 26 In Ref., 27 the Woods-Saxon basis was proposed as a 
reconciler between the harmonic oscillator basis and the coordinate space. The 
Woods-Saxon wave functions have more realistic asymptotic behavior at large r. 
One can also use a box boundary condition to discretize the continuum. It has been 
shown that to solve the spherical relativistic Hartree equations in a Woods-Saxon 
basis is almost equivalent to to do it in coordinate space. 27 The Woods-Saxon basis 
may be used in more complicated situations, e.g., for the description of exotic nuclei 
where both deformation and pairing have to be taken into account. 28 In this pro- 
ceeding paper we will present briefly the formalism, numerical procedure and some 
results of the deformed relativistic Hartree Bogoliubov model in a Woods-Saxon 
basis. Throughout the paper, this model will be labeled as the DRHBWS model. 

The paper is organized as follows. In Sec. [2 we give the formalism of the 
DRHBWS model. The numerical procedure and some preliminary results are pre- 
sented in Sec. [3] A summary is given in Sec. [4] 



2. Formalism 

The starting point of the relativistic Hartree theory is a Lagrangian density where 
nucleons are described as Dirac spinors which interact via the exchanges of several 
mesons (a, u>, and p) and the photon, 29-33 

C = V3i (ip - M) ^ + -d^a - U{a) - g a ^a^ 

- -fi^fi^ + irn 2 lo^ - g^ii/jipi 

- -R^R* 1 " + \^m 2 p p^ - gp^ipfijJi 

-ii^F^-e^^pM, (1) 

with the summation convention used, = "f^x^ = 7 M x M , M the nucleon mass, and 
m a , g a , m w , <? w , m pi g p masses and coupling constants of the respective mesons. The 
nonlinear self-coupling for the scalar mesons is crucial for a satisfactory description 
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of the surface properties 



34 



TTf \ 1 2 2 , 92 3 .93 4 
[/(cr) = -m a (T + y(7 J + — a , 



and field tensors for the vector mesons and the photon fields are defined as 
R»v = d^pv - dupfi - g P {p P x p v ), 



(2) 



(3) 



F„ 



ft A — ft A 

For the ground state of nuclei with time reversal symmetry, the nucleon spinors are 
the eigenvectors of the stationary Dirac equation, 

[a ■ p + V(r) + 0(M + S(r))} ^(r) = e^i(r), (4) 

and equations of motion for the mesons and the photon are 

(-A + d <7 U(a))a(r) = -g ( ,p s (r), 
(-A + ml) co°(r) = g uPv (r), 
(-A + mfj p°(r) = g p p 3 (r), 
-AA°(r) = epp(r), 

where u>° and A° are time-like components of the vector u and the photon fields and 
p° the 3-component of the time-like component of the iso-vector vector p meson. 
Equations (jl]) and (0) are coupled by the vector and scalar potentials, 



(5) 



V(r) = g uU °{r) + g P T 3 p°(r) + e^-A°{r), 
S(r) = g a cr(r), 



(6) 



and various densities 



Pv(r) = EiLiV'lWV'iW, 
Pc{r) = EiLiV'iW— ip-i>i(r)- 



(7) 



The pairing correlation is included via the Bogoliubov transformation and the 
relativistic Hartree-Bogoliubov equation reads, 35 

' h D (rap,ra'p') - A A(rap,r'a'p') \ / Uk(r'a'p') s 
—A*(rap,r'a'p') ~h D (rap,ra'p') + A J \ Vk(r'a'p') 



E 



d 3 r' 



( U k {rap)\ 



(8) 



\V k (rop) J ' 

where Ek is the quasiparticle energy, hp the Dirac Hamiltonian in and A the 
Fermi energy. Here p = 1,2 (or ±) is used to represent the particle-antiparticle 
degree of freedom. 

If in the pp channel, we use a zero range density dependent force, 



Vpi P2 p' 1 p' 2 i r l,r2;aia 2 a' 1 a 2 ) = jV 5{ri - r 2 ) [1 - Aa iv ■ a 22 >] P? x , • I^] > ( 9 ) 
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for axially deformed nuclei with spacial reflection symmetry, we may expand the 
potentials, S(r), V(r), and A(r), and various densities in terms of the Legendre 
polynomials, 

o \ i i r 

f(r)=^2h{r)P x (cos6), / A (r) = — — / d cos #/(r)P A (cost?) , A = 0,2,4,..-. 

A ^ 

(10) 

The quasi particle wave function is expanded in the Woods-Saxon basis 
{eiremj Vnnmirvp)} as, 



Uk(rap) = 

where the single particle energy ti Km and WcLVC function (pi^ui 

(rap) are obtained 

from solving the spherical Dirac equation with Woods-Saxon-like potentials, 36 

^ ' r \-F iK {r)Y} m {6, 4>,a) ) J 2' [ ' 

with (r) / r and i 7 ^ (r) /r the radial wave functions for the upper and lower compo- 
nents and Yj m (9, (f>, a) the spin spherical harmonics where k — (— l) J + i+1 / 2 (j + 1/2) 
and 1 = 1 + (— lp+'-V 2 . The states both in the Fermi sea and in the Dirac sea 
should be included in the basis for the completeness. In the Woods-Saxon basis, for 
each m-block the RHB equation (jHJ) turns out to be, 

^)("H(v)< 

where 



w=KSJ.v = fifty, (14) 



^=«U))' 2) =HsU))' (15) 

B= f A (™)_A c=('-A ( ~ ) =A (m) ~ ^ (16) 

\ (i K )(i'K')J ' V (m)(i'K') (i'«')(tre)/ ' 

The derivation of the matrix elements will be given in a detailed paper. 37 



3. Numerical procedure and results 

There are several parameters which have to be introduced for numerical reasons, 
e.g., the mesh size Ar, the box size R maK , and the cut off parameter on A in the 
expansion (|10[) . A ma x- Instead of a cut off on the radial quantum number n in the 
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Fig. 1. The total binding energy E (the lower panel), the proton pairing energy Sp air (the middle 
panel), and the rms radius R (the upper panel) of 20 Mg versus the cutoff energy E^~ ut in the Woods- 
Saxon basis for the DRHBWS model (solid circles). The spherical RCHB results (dashed lines) 
are also included for comparison. 



expansion (jTTJ) , an energy cutoff E^ ut is introduced for positive energy states in 
the Woods-Saxon basis and in each ft-block, the number of negative energy states 
in the Dirac sea is the same as that of positive energy states in the Fermi sea. 
We have investigated the dependence of our results on some of these parameters 
in spherical and deformed relativistic Hartree models. 27 ' 28 It's found that a box of 
the size i? max = 4r A 1 / 3 with r = 1.2 fm, a step size Ar = 0.1 fm, A max = 4, 
and E+ Ut = 100 MeV give relative deviations of the binding energy, the rms radius, 
and the quadrupole moment from the standard ones smaller than 0.1 % for light 
nuclei. In order to reduce the computational time, a small cutoff A max = 3 is used 
which would not introduce sizable errors. The parameter set NL3 38 is used for the 
Lagrangian density. 

There are two parameters (po and Vq) in the phenomenological pairing force ([9]). 
Since we are using a zero range force, a cutoff S^ut mus t be introduced to define 
the pairing window. We take the empirical value 0.152 fm~ 3 for the saturation 
density po- The pairing strength Vq — 374 MeV fm 3 and the cutoff E^' = 60 
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MeV reproduce the proton pairing energy in the spherical nucleus 20 Mg from the 
spherical relativistic Hartree Bogoliubov theory in a harmonic oscillator basis in 
which the Gogny force is used in the pp channel. 

As the first application of the DRHBWS method, we study a spherical nucleus 
20 Mg so that a comparison can be made between the DRHBWS and the spheri- 
cal RCHB 4,6 results. The ground state properties of 20 Mg are calculated with the 
DRHBWS and the RCHB codes. The total binding energy E, the proton pairing 
energy -Ep air , and the rms radius R of 20 Mg are plotted versus E* ut in the lowest, 
middle and top panels of Fig.Q] respectively. The spherical RCHB results, regarded 
as exact ones, are shown as horizontal dashed lines. When the basis size increases, 
E, -Spa;,., and R all converge to the corresponding exact values. In practical calcu- 
lations, one can choose -E^t according to the balance between the desired accuracy 
and the computational cost. For light nuclei, one can safely use E* ut = 100 MeV 
which results in accuracies in the total binding energy and the proton pairing energy 
of about a hundred keV and in the rms radius of around 0.002 fm. 

We are investigating some light deformed nuclei close to the neutron drip line 
by using the DRHBWS model. The results will be presented elsewhere. 37 

4. Summary 

In order to give a proper description of exotic nuclei, particularly deformed ones 
with large spatial extension, the deformed relativistic Hartree-Bogoliubov (DRHB) 
equations are solved in a Woods-Saxon basis in which the radial wave functions 
have proper asymptotic behaviors at large distance from the nuclear center which 
is crucial for the formation of halo. 

In this contribution, the formalism and the numerical procedure of the DRHB 
model in a Woods-Saxon basis are briefly presented. Some preliminary results, 
namely, the results for a spherical nucleus 20 Mg is also given and compared with 
the spherical relativistic Hartree-Bogoliubov model. 
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